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Superconductive quantum circuits (SQCs) comprise quantized energy levels that may be coupled 
via microwave electromagnetic fields. Described in this way, one may draw a close analogy to 
atoms with internal (electronic) levels coupled by laser light fields. In this Letter, we present 
a superconductive analog to electromagnetically induced transparency (S-EIT) that utilizes SQC 
designs of present day experimental consideration. We discuss how S-EIT can be used to establish 
macroscopic coherence in such systems and, thereby, utilized as a sensitive probe of decoherence. 
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Superconductive quantum circuits (SQCs) comprising 
mesoscopic Josephson junctions can exhibit quantum co- 
herence amongst their macroscopically large degrees of 
freedom 0. They exhibit quantized flux and/or charge 
states depending on their fabrication parameters, and the 
resultant quantized energy levels are analogous to the 
quantized internal levels of an atom. Spectroscopy, Rabi 
oscillation, and Ramsey interferometry experiments have 
demonstrated that SQCs behave as "artificial atoms" un- 
der carefully controlled conditions dH&SEHIllI|. 
This Letter extends the SQC-atom analogy to another 
quantum optical effect associated with atoms: electro- 
magnetically induced transparency (EIT) 0, [n|. We 
propose the demonstration of microwave transparency 
using a superconductive analog to EIT (denoted S-EIT) 
in a superconductive circuit exhibiting two meta-stable 
states (e.g., a qubit) and a third, shorter- lived state (e.g., 
the readout state). We show that driving coherent mi- 
crowave transitions between the qubit states and the 
readout state is a demonstration of S-EIT. We further 
propose a means to use S-EIT to experimentally probe 
the qubit decoherence rate in a sensitive manner. The 
philosophy is similar to that in Ref. IT^ . where it was pro- 
posed to use EIT to measure phase diffusion in atomic 
Bose-Einstein condensates. 

The three-level A system illustrated in Fig. HJi is a 
standard energy level structure utilized in EIT 0, ^| . 
It comprises two meta-stable states |1) and |2), each of 
which may be coupled to a third excited state |3). In 
atoms, the meta-stable states are typically hyperfine or 
Zeeman levels, while state |3) is an excited electronic 
state that may spontaneously decay at a relatively fast 
rate In an atomic EIT scheme, a resonant "probe" 
laser couples the |1) <-> |3) transition, and a "control" 
laser couples the |2) «-> |3) transition. The transition 
coupling strengths are characterized by their Rabi fre- 



quencies J7j3 = — dj3 • Ej3 for j = 1,2 respectively, 
where dj?, are the dipole matrix elements and Ej3 are the 
slowly varying envelopes of the electric fields. For par- 
ticular Rabi frequencies , the probe and control fields 
are effectively decoupled from the atoms by a destruc- 
tive quantum interference between the states of the two 
driven transitions. The result is probe and control field 
transparency In more recent experiments, ultra- 

slow light propagation due to EIT-based refractive index 
modifications in atomic clouds have also been demon- 
strated mom. 

SQCs have also been demonstrated to exhibit A-like 
energy level structures 00, 0,0- One example is the 
persistent-current (PC) qubit, a superconductive loop in- 
terrupted by two Josephson junctions of equal size and a 
third junction scaled smaller in area by the factor a < 1 
(Fig. [TJ)) |3 E3- Its dynamics are described by the 
Hamiltonian 

+ Ej [2 + a — 2 cos tp p cos <p m — a cos(27r/ + 2<p m )] , 

(1) 

in which C is the capacitance of the larger junctions, 
= ((f! ± (p2)/2, (fi is the gauge-invariant phase 
across the larger junctions i = {1,2}, Ej is the Josephson 
coupling energy, and / is the magnetic flux through the 
loop in units of the flux quantum <J>o- The second term 
in Eq. defines the magnetic-flux-dependent qubit po- 
tential landscape. For flux biases near one-half of a flux 
quantum, / « 1/2, the potential may be approximated 
by a double-well potential, with each well corresponding 
to a distinct, stable, classical state of the electric cur- 
rent, i.e., left or right circulation through the loop. In 
turn, each current state has a net magnetization of oppo- 
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FIG. 1: (a) Energy level diagram of a three-level A sys- 
tem. EIT can occur in atoms possessing two long-lived states 
|1), |2), each of which is coupled via resonant laser light fields 
to a radiatively decaying state |3). State |3) can feed back 
into |1), |2) and/or decay into levels outside the A configura- 
tion, (b) Circuit schematic of the PC qubit and its readout 
SQUID, (c) Schematic energy level diagram for a three-level 
superconducting quantum circuit. For our parameters we cal- 
culate U2 - wi = (2ir) 36 GHz and 1^3 — L02 = (2-7r) 32 GHz. 
The simulated matrix elements are (p|sin(27r/ + 2(/> m )\q) for 
(p,q) = (1,2), (2,3), and (1,3) are, respectively, 0.0704, -0.125, 
0.0158. 



site direction that is measurable using a dc SQUID [19j . 
As a quantum object, the potential wells exhibit quan- 
tized energy levels corresponding to the quantum states 
of the macroscopic circulating current 0, 0|. These 
levels may be coupled using microwave radiation 0, Il7| , 
and their quantum coherence has been experimentally 
demonstrated @- Note that for this system, the terms 
"population" and "occupation probability" are used syn- 
onymously. 

Tuning the flux bias away from / = 1/2 results in the 
asymmetric double-well potential illustrated in Fig. QJ;. 
The three states in the left well constitute the super- 
conductive analog to the atomic A system. States |1) 
and |2) are "meta-stable states," with a tunneling and 
coherence time much longer than the excited "readout" 
state 1 3). State |3) has weakly-coupled intra- well transi- 
tions, but has a strong inter-well transition when tuned 
on resonance with state |4) |16l 1 171 Hg| . Using tight- 
binding models with experimental PC qubit parame- 
ters [13, EE Ell at a flux bias f=0.5041, we estimate the 
tunneling times from states |1), |2), and |3) to the right 
well are 1/Ti w 1 ms, l/r 2 « 1 /is, and I/T3 w 1 ns 
respectively. Thus, a particle reaching state |3) will tend 



to tunnel quickly to state |4), and this event results in 
a switching of the circulating current that may be de- 
tected using a fast-measurement scheme. Alternatively, 
one may detune states |3) and |4), and then apply a res- 
onant 7r-pulse to transfer the population from state |3) 
to |4); since the states are now off-resonance, the relax- 
ation rate back to the left well is reduced and a slower 
detection scheme may be used. We note that a single- 
junction qubit shares this property, since the right well 
in Fig. nj; is effectively replaced by a quasi-continuum of 
states, and transitions out of the left well will not return. 

Transitions between the quantized levels are driven by 
resonant microwave- frequency magnetic fields. Assuming 
the Rabi frequencies f2,j to be much smaller than all level 
spacings \u>ki\ = \<-<Jk — the system-field interaction 
may be written within the rotating wave approximation 
(RWA)[2ll, 
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in which the decay from state |3) is treated phenomeno- 
logically as a non-Hermitian matrix element [2H |2^ . 
For small microwave perturbations of amplitude /a, 
the associated Rabi frequencies are given by Q pq = 
/a(p| sin (27r/ + 2<j) m )\q); numerical simulations of the 
matrix elements using PC qubit parameters are con- 
sistent with recent experimental results (see caption 
Fig. ^) 0, IrH Il8j . In general, all three intra- well tran- 
sitions are allowed in SQCs. For example, consider states 
|1) and 1 2) to be a qubit that is prepared in an arbitrary 
superposition state |\&) = ci|l) + ca|2) by temporarily 
driving the f2i2 transition. Then, by applying a probe 
field, the population of state |1) may be read out through 
a transition to state |3) followed by a rapid escape to the 
right well In this case, both the preparation and 
readout transitions were allowed and absorptive. 

One may achieve S-EIT in a superconductive A system 
that is prepared in state \^>) = ci|l) + c 2 |2), by simulta- 
neously and solely applying the microwave fields O13 and 
S7 2 3 such that 
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Under this condition (with SI12 = 0), the state \*&) is an 
eigenstate of tt^ r in Eq. with eigenvalue zero, and 
the SQC becomes transparent to the microwave fields. As 
in conventional EIT, the amplitudes for the two absorp- 
tion transitions into |3) have equal and opposite probabil- 
ity amplitudes, leading to a destructive quantum interfer- 
ence. Thus, in the absence of decoherence, preparing the 
qubit in state |\&) with an ideal preparation (^12) field 
and subsequently applying ideal probe (^13) and control 
(O23) microwave fields which satisfy Eq. © would result 
in no population loss through the readout state |3). In 
this way, S-EIT would confirm, without disturbing the 
system, that we had indeed prepared the qubit in the 
desired state. 
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In a practical SQC, there will be decoherence of the 
state and this must be measured, characterized, and 
minimized for quantum information applications. S-EIT 
is one such sensitive decoherence probe, since deviations 
in the amplitude and/or relative phase of the complex 
coefficients Ci from the condition established in Eq. @ 
result in a small probability |(cif2i3 + 02^23) /^| 2 of the 
SQC being driven into the readout state |3) on a time 
scale ~ r 3 /f2 2 . In general, there are two categories of 
decoherence: loss and dephasing. Loss refers to popu- 
lation losses from the metastable states |1), |2), and it is 
present in an SQC due to, for example, the finite loss rate 
of level |2>, T 2 - l/fxs = (2vr) 0.2 MHz. Dephasing refers 
to interactions of the SQC with other degrees of the free- 
dom in the system that cause the relative phase between 
ci and C2 to diffuse. The incorporation of dephasing is 
facilitated by the use of a density matrix formalism. 

We describe the system with a 3 x 3 density matrix 
with diagonal elements pa describing the populations, 
and pij, i ^ j describing the coherences between levels. 
In the presence of the EIT fields SI13 and ^23 with no 
direct coupling (Q12 = 0), the Bloch equations govern 
the evolution of the density matrix |2jJ : 



(a) 



% 1 

Pu = -Tipn - 2^*3P3i + 2^isPi3, 

f>22 = —F 2 p22 — 2^23^32 + 2^23/?23j 
P33 = -r 3 p 33 + |n* 3/ 31 - -fil3/3l3 

+ 2^23^32 - 2^23P23, 
Pl2 = -712P12 - 2^13^32 + 2^23Pl3, 

1 1 

Pl3 = -713P13 + 2^13(Pll - P33) + 2^23^12, 
P23 = -723P23 + 2^23(^22 - P33) + 2^13^21- 
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The remaining three elements' equations are determined 
by p*j = pji. The decoherence rates 7^ = (I\ + Tj)/2 + 
^(deph) j nc j u( j e both i oss anc [ dephasing contributions. 
We concentrate on the regime in which the readout state 
escape rate T 3 = 1 ns -1 = (2ir) 130 MHz dominates 
all other loss and dephasing rates, thus 713 » 723 ~ 
T3/2. Furthermore, we ignore the meta-stable state 
losses ri,r 2 relative to the dephasing 7|2° ph ^ an d se t 
712 ~ 7i2 Cph ' ) - Theoretical estimates of dephasing rates, 
such as , Yi^ eph \ in multi-level systems were recently ob- 
tained in Ref . |23j. 

We illustrate an S-EIT decoherence probe example by 
applying EIT fields fi 13 = ft 23 = (2tt) 150 MHz to the 
dark state \SS?) = (|1) — |2))/v / 2 with a dephasing rate 
712 = (2ir) 5 MHz and numerically integrating Eqs. J3J- 
. The dephasing 7 12 causes a small population transfer 
to the excited state P33 (Fig. EJl) . The excited state pop- 
ulation initially exhibits a rapid rise (see inset Fig. [2^,) 
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FIG. 2: (a) The population (occupation probability) P33 as a 
function of time with EIT fields Q 13 = n 23 = (2k) 150 MHz 
applied to an initial dark state pu = P22 = 0.5 and P12 = 
—0.5, and with a dephasing rate 712 = (2k) 5 MHz. The 
inset shows the same curve zoomed in on the early times. 
Under conditions of good EIT, we see a rapid initial rise to 
some plateau, followed by a much slower decay, (b) The 
total population P(t) remaining in the system versus time 
for the same simulation (solid curve). For comparison, the 
dashed curve shows the population for the out of phase case 
P11 = P22 = 0.5 and pi 2 = 0.5 discussed in the text. 



with transitory oscillations, reaching its maximum value 
p^ ay ^ within about 4 ns. This is followed by a smooth 
decay with a 1/e time of about 80 ns. The solid curve in 
Fig. (2t> traces the total population P = pn + p 2 2 + P33 
remaining in the system as a function of time. When 

the excited state maximum p^ 1 "^ is reached, the total 
remaining population is P(4 ns) = 0.973. In contrast, 
the dashed line in Fig. \2)p illustrates the rapid popula- 
tion loss expected when the same fields are applied to 
the state |*) = (|1) + |2))/V2 [ir out of phase with the 
dark state in Eq. (J3J]. In the absence of S-EIT quan- 
tum interference, the entire population is lost on a time 
scale T^/fl 2 ~ 1 ns. The general behavior presented in 
Figs. [2^ and[2b is observed over a wide parameter regime 
of experimental interest. 

We now use Eqs. to show how measuring the 

slow population loss in S-EIT can be used to extract the 
decoherence rate 712. The elements P33, P13, and P23 
in Eqs. El El and El are damped at a rapid rate ~ T3, 
allowing their adiabatic elimination |22j, Il24| ; we solve for 
their quasi-steady state values by setting P33 = P13 = 
P23 = 0. This approximation is accurate once initial 
transients have passed and the plateau value p^^ has 
been reached. Using these results in Eq. J7J) yields an 
equation for P12 with a strong damping term fl 2 /T3, and 
it too can be solved for its quasi-steady state value. In 
the limit 7i 2 r 3 /0 2 «lwe get [U 
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(Pll(t)+P22(t)). (10) 



The ratio 27i2F3/i7 2 represents the small fractional de- 
viation of P12 from its dark state value. There is a com- 
petition between the "preparation rate" VL 2 /T3 (which 
constantly acts to drive the system into the dark state) 
and the decoherence rate 712 (which attempts to drive it 
back out). 
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FIG. 3: (a) The maximum plateau value P3™ ax ' for different 
712 (circles). The solid curve shows the prediction (b) 
The remaining population P — p\\ +P22 + pzz at the time the 
plateau is reached for the cases in (a). 



We now consider the conditions under which one can 
use S-EIT to estimate the decoherence rate with little 
population loss in the system. Eqs. JSJ, and i|l(J|) 

reveal that deviations from the dark state cause popu- 
lation loss through |3) at a rate R = 27i2(£li 3 f2 2 ! 3 /f2 4 ); 
this ultimately leads to the exponential decay of P seen in 
Fig. [2t>. By assumption, population escapes the system 
only through the decay term — p 33 r 3 in Eq. yielding 

P33 = {R/^?>)P- At the early time when p 3 ™ has been 
reached, the population P remains close to unity, and the 
excited state population reaches 
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(11) 



The time T ss to reach p 3 ™ ax ' ) is generally the smaller of 
the preparation time ~ r 3 /f2 2 and the inverse of the de- 
cay rate I/T3. At this time, the total population loss will 

be ~ t ssP { 3 T x) ~ (2n 2 13 nl 3 /n 4 )M a x( 112 /T 3 , 112 r 3 /n 2 ). 

So long as the loss during this initial transient time is 
small, the population will follow a simple exponential 
decay P(t) — exp(— p 3 ™ ax ^r 3 i), and the dephasing rate 
712 can be easily extracted. To keep this loss small, we 
require both ratios in the Max(- • • ) argument to be small 
(fl ^> V27i 2 r 3 and T 3 ^> 712) in order to use this ap- 
proach to accurately estimate the decoherence rate while 
causing little loss from the system. Since f2 is experi- 
mentally controllable, it can be chosen to satisfy the first 



constraint. If T 3 is comparable or smaller than 712, then 
S-EIT remains a decoherence probe, although the strong 
damping assumption leading to Eq. (|lf)|> no longer holds 
and so the analysis is different. 

We have performed a series of numerical simulations, 
varying 712 to test the validity of the above approach. 
The results are presented in Fig. [3J Fig. [SJi indicates 
Ps® 3 *^ versus 712 and compares the results with the 
analytic estimate [Eq. (|llfl ]. The agreement is quite 
good for 712 < (2ir) 4 MHz, which corresponds to the 
27i2r 3 /f2 2 < 0.056. Higher dephasing rates compete 
more with the preparation rate, making the adiabatic 
elimination approach less valid; this leads to deviations 
from our analytic prediction [Eq. (fill) ]- In such cases, one 

observes a significant loss by the time p 3 ™ ax ' is reached, 
as illustrated in Fig. |2t> 

The RWA invoked in the calculation ignores far off- 
resonance couplings induced by other applied fields at 
other transitions (e.g., in our case, the 5123 field drives 
|1) «-> 1 2) at A ~ 4 GHz off resonance). We have per- 
formed calculations including all such off-resonant cou- 
plings. Generally, they lead to small shifts of the en- 
ergy levels (analogous to A.C. Stark shifts) and loss 
rates. These losses scale as T 3 (fl 2 /A 2 ) and fl 4 /T 3 A, 
putting a limit on the field strengths fl 2 which can be 
used for a given A. For the parameters considered here, 
we found shifts of (2ir) 8 MHz and loss rates totalling 
(2?r) 19 kHz; this should not effect measurements in the 
regime 712 > (2tt) 19 kHz j^. 

We have proposed using the superconductive analog 
to EIT (S-EIT) to demonstrate macroscopic quantum 
interference in superconductive quantum circuits. S- 
EIT provides an accurate and sensitive means to probe 
the accuracy and phase coherence of qubit preparation, 
and we have calculated analytic expressions for the field 
strengths required for this purpose. 
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